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04 ' We dedicate the paper to the memory of our dear friend and colleague Juha Heinonen (1960-2007). 
Ft . ■ Abstract. Let / : ^ be a mapping of finite distortion, where Q C M^. Assume 



(N 



that the distortion function K{x,f) satisfies e^^'-'^^^ e Lf^^{fl) for some p > 0. We 



lo 

estabUsh optimal regularity and area distortion estimates for /. In particular, we prove 
that l-D/P log'^~^(e + \Df\) G Ll^^{fl) for every (3 < p. This answers positively, in 
dimension n — 2^ the well known conjectures of Iwaniec-Sbordone [18. Conjecture 1.1] 
\^ \ and of Iwaniec-Koskela-Martin fT6l Conjecture 7.1]. 

' 1. Introduction 

We say that a mapping / = defined in a domain i7 C is a mapping of finite 

CN ■ distortion if 

^ i (i) / e 

(N . (ii) J = u^Vy - UyV^ e Lj^^i^), and 



(iii) tliere is a measurable function K{z) > 1, finite almost everywhere, such that 
\Df{z)\'^ < K[z) J{z, f) almost everywhere in Q. 



^ ' The smallest such function is denoted by K{z, f) and is called the distortion function of 
> > _ 

In two dimensions the mappings of finite distortion are intimately related to elliptic 
^ . PDE's. For equations with non-smooth coefficients the first condition, requiring that / has 

locally integrable distributional partial derivatives, is the smallest degree of smoothness 
where one can begin to discuss what it means to be a (weak) solution to such an equation. 

The second condition is a (weak) regularity property which is automatically satisfied 
by all homeomorphisms / in the Sobolev class Wl^liVt). The last condition, that the 
distortion function 1 < K{z^ f) < 00, merely asks that the pointwise Jacobian J{z, /) > 
almost everywhere and that the gradient Df{z) vanishes at those points z where J{z, f) = 
0. This is a minimal requirement for a mapping to carry geometric information. 
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In two dimensions the distortion inequality 
notation as 



df 



dz 



< k{z) 



df 



dz 



where k{z 



can be reformulated with the complex 

K{z) - 1 



K{z) 



1 



< 1, 



using the identifications Df = \ fz\ + \fz\ for the operator norm and J{z, /) = l/zP — 
for the Jacobian. Yet another equivalent formulation is with the Beltrami equation, 



(2) 



df , .df K{z) - 1 ^ ^ 



where the coefficient /i can be defined simply via the differential equation ([2]). Then at a 
given point z we have K{z, /) < oo if and only if \[ji{z)\ < 1. 

The classical theory of quasiconformal mappings considers precisely the mappings of 
bounded distortion, that is mappings with K{z, f) G or, equivalently, mappings with 
complex dilation \^{z)\ < k for some constant k < 1. The classical measurable Riemann 
mapping theorem tells that if we have a coefficient < k < 1 for 2; G C, then 

the equation ([2]) has a homeomorphic solution and all other solutions are obtained by 
post-composing with holomorphic functions, see [8], [1], or [H Chapter 5]. 

It is natural to ask how much the condition K{z, f) G can be relaxed in order 
to obtain a useful theory. In the groundbreaking work [10] David generalized the mea- 
surable Riemann mapping theorem to maps of exponentially integrable distortion, where 
in general we have ||/i||oo = 1- Later, Iwaniec, Koskela, and Martin [16] generalized 
the corresponding regularity theory to higher dimensions, and for the two dimensions a 
systematic modern development was established by Iwaniec and Martin as part of their 
monograph [T3]. During the last 10 years there has been an intensive study of finite distor- 
tion mappings, motivated also by the fact that these maps have natural applications e.g. 
to non-uniformly elliptic equations and elasticity theory. We refer the reader for instance 
to [ll],[15], [9], [23], [16] or [12] and the references therein, for the basic literature on the 
subject. 

As a typical example of a mapping with finite but unbounded distortion consider the 
following exampl^. Let 



(3) 



A/ 



log I e + 

\z\ 



n -p/2 



log log e + —7 

\z\ 



n-1/2 



for \z\ < 1, 



and for |z| > 1 set gp{z) 
that 



cq z. For each value of the parameter p > one easily computes 

while e^^^'^p^ is not locally integrable to any power s > p. Furthermore, we have the 
fundamental regularity Qp G W^;of(C) when p > 1. However, for p = 1 the mapping 

From example ([3]) we see that the integrability of e^*^^^ is in general not sufficient for 

1 2 

the Wi^^ -regularity. Instead we have a refined logarithmic scale of regularity around the 



""^The similar formula (20.75) in |4] has a slight misprint. 
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space Wl^^. Namely, for each p > 0, 

|D^p|2log^-i|D/|GLL(C), 0<P<p, 

while the inclusion fails for (3 = p. 

The exponential integrability of the distortion function in this theory is no coincidence. 
Indeed, to ascertain the continuity of the mapping / of finite distortion, unless some addi- 
tional regularity is required, the exponential integrability of K suffices but this condition 
has only little room for relaxation, see Theorem 11.2.1 and [4], Section 20.5. Hence 
the mappings of exponentially integrable distortion provide the most natural class where 
to look for a viable general theory. 

Within this class of mappings one of the fundamental questions is their (optimal) reg- 
ularity. That is, given a homeomorphism with a'priori only / G Wl^^{fl), we assume 
that 

^K{zj) ^ i:P^(fi)^ where p > 0, 

and then seek for the optimal improvement on regularity this condition brings, i.e. find 
in terms of p the best possible Sobolev class where the mapping / belongs to. 

The example ([3]) suggests natural candidates for the extremal behaviour. Indeed, in 
their paper Iwaniec and Sbordone [TH] formulated a precise general conjecture on the 
regularity of mappings with exponentially integrable distortion. The main purpose of this 
paper is to give a proof to their conjecture, and establish the following theorem. 

Theorem 1.1. Let Q G M."^ be a domain. Suppose that the distortion function K{z, f) of 
a mapping of finite distortion f G Wl^^{Q) satisfies 

(4) e^(^'^) G Ll^{n) for some p > 0. 
Then we have for every < P < p, 

(5) \Df\' \og^-\e + \Df\) G L^n) and 

(6) J(^,/)log^(e + J(;.,/))GLL(^^). 

Moreover, for every p > there are examples that satisfy yet fail and for 
P = p. 

Of particular interest is the question of when the derivatives of the mappings are L^- 
integrable. As a special case of our theorem one has 

Corollary 1.2. Let Q G be a domain and f G Wl^^{Q) a mapping of finite distortion. 
Suppose that 

(7) e-^(^'^) G LL(^) for some p > 1. 
ThenfeWl^^in). 

The family of mappings shows that W^^^ -regularity fails in general at p = 1. 
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Previously |T6], [18] it was known that there exists some constant po > 1 so that e^^^'^^ G 
-^foc with p > po imphes / G W^^^. Our result identifies the precise regularity borderline 
Po = 1. Similarly [H] (see [17], [12] for ra-dimensional results) the conclusions ([5]) and 
([6]) were only known to hold under the stronger assumption (3 < Cip, where Ci G (0, 1) is 
some unspecified constant. 

Our proof of optimal regularity employs the approach of David [10] . It turns out that 
it is practically impossible to describe our proof without delving in depth into David's 
mapping theorem. Hence, to make the exposition as clear as possible we will include a 
straightforward proof of this result, and present this in the modern framework developed 
by Iwaniec and Martin [14J: 

Theorem 1.3. [David, Iwaniec-Martin] Suppose the distortion function K = K{z) is 
such that 

e^ G L^(d) for some p > 0. 

Assume also that ^{z) = for |2| > 1. Then the Beltrami equation fz{z) = ^{z) fz{z) 
admits a unique principal solution f for which 

/G<f(C), Q{t)=t'\og-\e + t). 

The principal solution is a homeomorphism. Moreover, every other Wl^^ -solution h to 
this Beltrami equation in a domain Q G C admits the factorization 

/i = o /, 

where (j) is a holomorphic function in the domain /(f^). 

By a principal solution we mean a homeomorphism / G W^^^ (C) that satisfies 
(8) = , + ^ + ^ + ... 



z z"^ 



outside some compact set. 



There is a simple explanation why it is convenient to develop the theory of mappings 
of exponential distortion using the space Wl^^ as a starting point. Namely, suppose 
/ G Wi^^ (c) is an orientation preserving homeomorphism, that is J{z, f) > 0, whose 
distortion function K{z) satisfies e^ G L^^^{c) for some < p < oo. The Jacobian of any 
planar Sobolev homeomorphism is locally integrable. Hence we may use the elementary 
inequality 

ab < alog(l + a) + e'' - 1 

to find that 

log(e+|D/|2) - log(e + J) - p log(e + J) - p^ > 
for all p > 0. Thus 

for any bounded domain. This shows that / automatically belongs to the Orlicz-Sobolev 
class Wlf{C). 
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Conversely in \IA\, Theorem 7.2.1 it is shown that for an orientation preserving mapping 

/ the log~^ L-integrabihty of the differential Df implies that J{z, f) G Lj^^i^il). This 
slight gain in the regularity of the Jacobian determinant is precisely what makes the 
space Wl^^ the natural framework for the theory of mappings of exponentially integrable 
distortion. 

The structure of this paper is the following. First in Section [2] we recall the basic 
modulus of continuity estimates for mappings of finite distortion. Section [3] presents the 
crucial sharp estimates for the rate of decay of the Neumann series of the solution to 
the Beltrami equation, in the case of an exponentially integrable distortion. Theorem 
11.31 is then proven in Section HI where our presentation is self-contained, apart from the 
generalized uniqueness and Stoilow factorization theorem due to Iwaniec and Martin (see 
Theorem 14.21 below) . The proof of our main result. Theorem 11.11 is the content of Section 
O In that section we also obtain optimal area distortion estimates. The optimality of 
our results on the rate of convergence of the Neumann series is considered in Section 
El Finally, Section [7] gives an application of the optimal regularity results to degenerate 
elliptic equations. 



There is a particularly elegant geometric approach to obtaining modulus of continuity 
estimates for functions that are monotone. We will not consider the notion of monotonicity 
in its full generality, but just say that a continuous function u in a domain Q is monotone, 
if for each relatively compact subdomain Q' G Q we have 



that is, if u satisfies both the maximum and minimum principle. For instance, coordinates 
of open mappings are monotone. 

The idea of the next well known argument goes back to Gehring in his study of the 
Liouville theorem in space. 

Lemma 2.1. Let u G W'^''^{2n) he a continuous and monotone function. Then for all 
a, b & ID we have 



2. Modulus of Continuity. 



maxu = maxM and 
an' n' 



mmu = mmti, 

an' n' 



(10) 



n(a) -n(6)p < 



TT 




Proof. Consider the two concentric disks IS>{z,r) and I]){z,l), where z = 
r = ||a — Note that all disks ID^z, t) with t < Pq = max{l, 2r} lie in 3B. 
We have by monotonicity 



|(a + 6) and 




for almost all r < t < 1. By Cauchy-Schwartz 
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We integrate with respect to t, r < t < po, to get 



\u{a) — u{b)\'^ \og{po/r) < — / \Vu\'^dm < — j |VMp(im. 

2 Jn{z,po) 



The estimate (|T0|) is a quick consequence. □ 

In view of this result uniform bounds for the L^-derivatives become valuable. Within 
the (uniformly elliptic) Beltrami equation such bounds can be obtained for the inverse of 
the principal solution. Iwaniec and Sverak [19] were the first to make a systematic use of 
this phenomenon. 

Theorem 2.2. Let f G Wl^^{C) be the principal solution to fz{z) = fi{z)fz{z), z E C, 
where p is supported in the unit disk B with \\p\\oo < 1- Let g = f^^ be the inverse of f . 
Then 



[ i\9w\^ + \9..-l\^)dw<2 [ K{zJ)dz. 



Proof. Since the function g{z) — z has derivatives in Iv^(c), an integration by parts shows 
that the integral of its Jacobian vanishes. This gives 



gw-i\ - \gw\ = / Jiw,g) = o 



g^-l\^ = / |(?^p= / \g^\^< / \Dg\^ 

C Jc Jf{0) "'/(D) 



Therefore 



K{w,g)J{w,g)= / K{f{z),g)= / K{zJ) 

/(D) 

since quasiconformal mappings satisfy the usual rules of change of variables. □ 

The key fact in the following modulus of continuity estimate, as in Theorem 12.21 is 
that the bounds do not depend on the value of ||Ai||oo- It is obtained easily by combining 
Lemma [2. II with Theorem 12. 2[ 

Corollary 2.3. Let f and g = f^^ be as in Theorem \2.2\ and let R> 1. Then 



\a-b\ 

whenever a,b E B(i?). 

3. Decay of the Neumann Series. 

According to the classical measurable Riemann mapping theorem [Hj , [1] , given a com- 
pactly supported Beltrami coefficient with \fi{z)\ < k < 1 almost everywhere, the system 

always admits a homeomorphic solution / with the development ([8]). That is, / is a 
principal solution. 
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This uniformly elhptic equation is most conveniently solved by applying a Neumann- 
series argument using the Beurling operator S. This is defined by the principal value 
integral 

The Beurling operator acts unitarily on L^(c). Also, we have S{hz) = for every 
h e W^'HC). 



f = z + Ciuj), Cujiz):=- [ 



We then look for the solution in the form 

(z-t) 

where C is the Cauchy transform and uj = fz is to satisfy the identity 
uj{z) = fi{z)Suj{z) + fi{z) for almost every z E C. 
Such an u is easy to find when ||/i||oo < k < 1. We define 

(11) U := (1 — flS)~^fi = II + yU5/i + flSflSfl + fiSfiSfiSfl + ■ ■ ■ 

The series converges in L^(c) since the operator norm of the n:th iterate 

1 1 /i5/i>S ■■ -/iiS II L2(c)^L2(c) < fc", n G N. 

The converging series provides us with the solution /. However, showing that / is a 
homeomorphism takes some more effort. For this and other basic facts in the classical 
theory of quasiconformal maps we refer the reader to [T], [1] and f2T\ . 

It is natural to study how far the above method carries in the situation of degenerate 
equations, even if now ||/i||oo = 1 and the Neumann series of (I — fiS)~^ cannot converge 
in the operator norm. Hence the main question here is the rate of decay for the L^-norms 
of the terms in the series f fTT]) . The next result, building on the approach initiated by 
David \lQl, gives an optimal answer. 

Theorem 3.1. Suppose \^{z)\ < 1 almost everywhere, with ^{z) = for \z\ > 1. If the 
distortion function K(z) = satisfies 

e^'GLP(D), p>0, 
then we have for every < j3 < p, 



|(/i5)>r < Co(n + l)-^ nGN, 

where Co = Cp^p ■ e^^ with Cp^p depending only on (3 and p. 
Proof. Since 

K{z) + 1 = ^— ^, 

1 - \ii[z)\ 

with Chebychev's inequality we have the measure estimates 



|{2 e D : \^i{z)\ > 1 - i}| < e-^P' [ e*'(^+^) = Ce-'^* 
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for each t > 1. With this control on the Beltrami coefficient we will iteratively estimate 
the terms {fiS)"'fi of the Neumann series. For this purpose we ffist fix the parameter 
< (3 < p, and then for each n G N divide the unit disk into the "bad" points 

Br, = {zeB:\fx{z)\>l-^-^} 

and the "good" ones, i.e. the complements 

The above bounds on area read now 

where Ci = Jj, e^^ . 

Next, let us consider the terms ipn = of the Neumann series, obtained induc- 

tively by 

It is helpful to ffist look at the following auxiliary terms 

9n = Xg„ ^S{gn-i), go = /i, 

that is, in each iterative step we restrict /i to the corresponding good part of the disk. 
The terms gn are easy to estimate, 

2 



Thus 



where = a/vt (1 + /5/2)^^^. For the true terms ipn decompose 
This gives the norm bounds 

(12) \\^n-9n\\h< (l-^;^) ||^„_i-(7„_i||i. + i?(n), 

where 



i?(n) = |ixB„/i5(^„-i)lli. = / |(/i5)>|^ 

J Bn 

A main step of the proof is to bound this last integral term. A natural approach is to 
apply Holder's inequality and estimate the norms || (/i5)"/i||p for p > 2. This approach 
can be carried over by applying the spectral bounds on the operator fiS obtained in [5]. 
On the other hand, these bounds are based on the the area distortion results from [2]. 
To reveal more clearly the heart of the argument we will therefore directly use the area 
distortion. This can be done with a holomorphic representation as follows. 

Recall that the principal solution f = to the equation 

(13) Mz) = \fiiz)Uz) 



OPTIMAL REGULARITY FOR PLANAR MAPPINGS OF FINITE DISTORTION 9 

depends holomorphically on the parameter A G ID). This fact follows from the Neumann 
series representation for the derivative, 

(14) = A/i + x^fiSfi + ■■■ + ri^isy-^^i + ■ • • 

The series converges absolutely in L^(c) since ||A/i||oo < |A| < 1 and S is an L^-isometry. 

The L^(c)-valued holomorphic function A —» can as well be represented by the 
Cauchy integral. From this we obtain the following integral representation, 

(15) XE {f^STf^ = ^ I -T^ftxE dX, ECB, 



2m A«+2 

valid for any < p < 1. We are hence to estimate the norms 

II/.' XE\\h = I \m' < 1 Ji^' /') = l/'(^)l' 

and it is here that the need for a quasiconformal area distortion estimate arises. We have 
the uniform bounds 

(16) \f\E)\<7rM\E\'/'', 1^1 = ^^' M>1, 

see [6], Theorem 1.6 or [Ij, Theorem 13.1.4. 

If we now denote p = |A| = and combine f|T5|) with the above estimates, we end 
up with the result 



(17) 



|,,(,5)>||.<v^(|±i)"^£±i|^|V(-.. 



The estimate is valid for every M > 1 and for any Beltrami coefficient with < Xa 
almost everywhere. 

For later purposes we recall that the 9^-derivative, too, has a power series representa- 
tion, 

f^-l= Sf^ = XSfi + X'^SifiSfi) + ■■■ + A"5(/x5)"- V + • • • 
With an analogous argument we obtain 

(18) ||,,5(,5)>||,<V^(^)""^|i^l'"-'. 

an estimate which is similarly valid for every M > 1. 

Let us then return to estimating the original Neumann series. Unwinding the iteration 
in (fT2|) gives 

j=l k=j+l ^ j=l 

It remains to show that the remainder term R{n) decays exponentially. But here we may 
use the bound (fTTl) for the set E = Bn- We have |-B„| < Ci e"^"^/'^ and hence 

2n / » ^- , 1 \ 2n 



^ ' - \M-\) ' ' - ^ \M-\ 



_ 4n p 

6 ? , 



10 



K. ASTALA, J. GILL, S. ROHDE, AND E. SAKSMAN 



Given [3 < p we can choose M > 1 so that 

for some S > 0. With this choice R{n) < Ce-^''", n G N, where C = AM'^e'P j^eP^ . This 
completes the proof. □ 

According to the Theorem, the terms of the Neumann series decay with the rate 

||(yu5)"/i||2 < C/3r^-^/^ for any p<p. 

If the decay would be a little better, of the order ra"^, then for any p > 1 the series 
would be norm convergent in L^(c). In view of the example ([2]), this would immediately 
determine po = 1 as the critical exponent for the 1^^' ^-regularity. However, as we will see 
later in Section [6], the order of decay in Theorem 13.11 cannot be improved. Hence further 
means are required for optimal regularity, to be discussed in Section [51 

Nevertheless, at the exponent po = 1 there is an interesting interpretation of the above 
bounds, in terms of the vector-valued Hardy spaces. Namely, suppose /i is as in Theorem 
Oandp > (3>l. Then 

oo 

(19) 5^||(/.5)>f^^<oo. 

n=0 

On the other hand, as in ([T3D and dTl), consider A ^ holomorphic vector val- 

ued function D L'^{c). With the above bounds the function has square-summable 
Taylor coefficients {fiS)"'fi. In other words, it belongs to the vector- valued Hardy space 
i7^(D; L^(c)). In particular, since the target Banach space L^(c) has the Radon-Nikodym 
property, the radial limits 

lini^eL2(c), Ces\ 

exist in L^(c) for almost every G see e.g. [7]. Similarly for the derivative dzfr(;, and 
it is not difficult to see that the limiting derivatives are derivatives of a principal solution 
/(;, solving 

(20, t-C-wf- 

m 1 12 

Thus already here we obtain, for p > 1 and for almost every G § , the VF^J^ -regular 
solutions to (12(1 . 

For exponents /3 > 2, we have via the norm convergence (|T9l) the next auxiliary result, 
which the reader may recognize as a distortion of area for mappings in the exponential 
class. The result here is rather weak but is a necessary step towards the optimal measure 
distortion bounds which will be established in Section [5l 

Corollary 3.2. Suppose /i and < P < p are as in Theorem \3.1\ and 

oo 
n=0 
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If 13 > 2, then 

\\XEO^\\2+\\xESa^\\2<C\oi-^'\e + ^^, ECB, 
where C = CqC^ < oo with Cp depending only (3. 

Proof. Theorem 13.11 and its proof give us two different ways to estimate the terms. First, 

||xE(/i5)>||2 + \\xES{^lSr^lh < '^WSTi^h < 2Co{n + i)-^'\ 

Summing this up gives 

n=m+l 

Secondly, let us choose, say, M = 3 in (1171) and fll8l) . Then for every n G N we have the 
estimates ||x£;(/^'5)"'/x||2 + \\xe<S{ij^S)'^h\\2 < ■ 2" ji?!^/^. This, in turn, leads to 

m 

J2 \\XE{f^S)yh + WxESifiSr^ih < 6v^ • 2"+i|E|i/^ 

n=0 

Combining we arrive at 

WXE + WXE Sa,h < ^ rn'-'^' + 22 ■ 2'"|i?|V6. 

The bound holds for any integer m G N, but of course it is optimal when the two terms 
on the right hand side are roughly equal. We let 

— log(e + ^) ^ < 1 + log(e + 

and with this choice the terms on the right hand side of (!2T!) are both smaller than 
Cq Cp log^~^^^(e + j^), where the constant Cp depends only on the parameter (3 > 2. □ 



4. The Degenerate Measurable Riemann Mapping Theorem 

We will next give a simple proof of David's theorem, the generalization of the Measur- 
able Riemann Mapping Theorem to the setting of exponentially integrable distortion. We 
give this presentation not only for the reader's convenience, but also since a few bits and 
pieces of the proof will be needed in the main result, the optimal regularity Theorem II. 1[ 

We start with the case where the exponential distortion e^ G Lf^^ for some p > 2. The 
key idea here is that for such large p, the Neumann series (fTTj) converges in -L^(C), and 
this quickly yields homeomorphic ly^'^-regular solutions to the corresponding Beltrami 
equation ([2]). From Section [2] we already know the usefulness of such regularity, providing 
us with general equicontinuity properties. 

After the case of high exponential integrability, together with a few further consequences 
established, we then return to the general measurable Riemann mapping theorem and 
prove this towards the end of this section. 
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Theorem 4.1. If ji is a Beltrami coefficient such that 

, , Kiz) - 1 

where 

e LP(ID)) for some p > 2, 

then the Beltrami equation 

df , ,df , , 

-— = u[z)—— for almost every 2; G C 

admits a unique principal solution f G Wlj^iC). 

Proof. Consider the following good approximations for the Beltrami coefficient ^{z), 

r^-J-^^^) if|M^)|<i-^ 

^^^> l^m{z)-< M(£) otherwise, 

defined for m = 1, 2, . . .. Of course < 1 — ^ < 1 and we also have the bound 

K[z) + 1 

which is independent of m. Similarly the corresponding distortion functions satisfy 

A-,„(.) = < = K(z), m = 1. 2. . . . 

Note that each Km{z) G L°°(c). 

With the classical measurable Riemann mapping theorem, see [8], [1], [1] or [21], we 
have unique principal solutions : C — C to the Beltrami equation 

(n. = /^m(n., m = i,2,... 

Moreover, 

f^ = Z + C{Um), 

where ujm = {.f^)z satisfies the identity 

^m{.z) = fim{z){SuJm){z) + fim{z) almost cvcry z G C. 

In order for these approximate solutions to converge we need uniform L^-bounds for 
their derivatives; once these have been established the proof follows quickly. 
We use Theorem 13.11 This gives for any 2 < (3 < p the estimates 

II (/i„5)Vm II 2 < ken. 

The theorem also gives uniform bounds for the constant term Cm- Since the approximate 
solutions have distortion Km < K pointwise, we have 
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where Co is as in Theorem l3.1[ Hence if G L^(d) for some p > 2, we obtain L^-estimates 
for the sums (fTTj) . 



(22) ||(n.l|2<X]||(/^m5)Vn^l|2< V^oJ^i^ + l] 



k=0 k=0 

Hence also 

oo 

-(3/2 



Ii(n.-iii2 = ii(n.ii2< v^E(^ + i) 

k=0 

This easily yields that 

(23) / \Dr\^ KAttR"^ + A, R>1, 

Jo{R) 

where A < oo depends only on /5 > 2 and e^^ . 

Now Lemma 12.11 applies and shows that on compact subsets of the plane the sequence 
Ijmj j^g^g ^ uniform modulus of continuity. Changing to a subsequence, and invoking 
( |23|) again we may assume that f"^{z) f{z) locally uniformly with the derivatives Df"^ 
converging weakly in Lfg^{c). Their weak limit necessarily equals Df. Furthermore, since 
for any > 1 and ^ e C^{B{0, R)) 

[ v^(/r-/i/r)= / Hf'm - ^) fT < ^ V47ri?2 + A ^ 

Jc Jc 
by (l23|) . we see that / is a solution to the Beltrami equation 

fz=fJ'fz for almost every z E C. 

On the other hand, we can use Corollary 12.31 for the inverse mappings Qm = (Z™)"^ and 
see that g"^{z) g{z) uniformly on compact subsets of the plane, where g : C C is 
continuous. From g"^ {f"^{z)) = z we have g o f[z) = z. That is, / is a homeomorphism, 
hence the principal solution we were looking for. 

Lastly, for the uniqueness we refer to the following very general Stoilow factorization 
result from the monograph [Hj; see also [1], Section 20.4.8. The proof of Theorem 14.11 is 
now complete. □ 



Theorem 4.2. ([H], Theorem 11.5.1). Let 

(24) Qit) 



log(e + 1) 

and suppose we are given a homeomorphic solution f G Wl^^{Q) to the Beltrami equation 
(25) /, = /i(z)/„ z G fi, 

where |/i(-2)| < 1 almost everywhere. Then every other solution h G Wl^^{Q) to (1251) takes 
the form 

h{z) = (Pifiz)), z G n, 

where (p '■ /(^) ^ C is holomorphic. 
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In Theorem 14.11 we constructed the mapping / through a hmiting process rather than 
directly by the Neumann series. Nevertheless, the representations 

oo oo 

(26) /, = Y,{fiS)y, f.-l = Y. 

n=0 n=0 

are still valid. Namely, by Theorem 13. li the sums are absolutely convergent in L^(c). Since 
the approximate coefficients fim — > in L°°(c), the terms {nmSYum {f^<S)''f^ weakly in 
L^(c). As the derivatives of /™ are representable by the corresponding Neumann series, 
the identities fl2Bl) follow with the help of the uniform bounds f l22l) . 

In studying the general exponentially integrable distortion we need to know that the 
mappings of Theorem 14. II preserve sets of Lebesgue measure zero. This follows easily from 
the above. 



Corollary 4.3. If fi and the principal solution f G W^jof (C) are as in Theorem 4jJ_. then 
f has the properties Af and J\f~^, i.e. for any measurable set E G C, 

|/(E)|=0 ^ \E\=0. 

Proof. Any homeomorphism with the Wl^^-regulahtj preserves Lebesgue null sets, see 
[22j . or j4j for an elementary proof. We have shown that / G Wi^^ (C), and since K{-, f) G 



L^(d), Theorem 12.21 verifies that the inverse map g = h ^ belongs to VT/of (C). The claim 
follows. □ 

Concerning the existence of solutions, the situation is rather different in the general 
case where the integrability exponent of is small. Recall that by example ([3]) the 
corresponding principal solutions need not be in VF^of (C). Instead, inequality ([9]) leads us 
to look for principal solutions / in Wl^^{c), where Q{t) is given in fl2^ . The idea of the 
following proof is then to reduce the argument to Theorem 14. II via a suitable factorization. 



Proof of Theorem M.'A Given a Beltrami coefficient /i = ^{z) with the distortion function 



^i.^) = Trrrrl) assume that G U'ip) for some < p < 2. It is very suggestive to 



write 

3 pK{z) 



K{z) 



p 3 

Then Ki{z) = pK{z)/3 satisfies the hypotheses of Theorem 14.11 and the constant factor 
K2{z) = 3/p > 1 could be represented as the distortion of a quasiconformal mapping. 
There is one problem here in that at points where K{z) is already finite and perhaps small, 
the factor Ki{z) might be less than 1 and so cannot be a distortion function. We will get 
around this point by constructing the related Beltrami coefficients using the hyperbolic 
geometry. 

Let M > 1. For each z ^ D choose a point u = v{z) on the radial segment determined 
by so that 

1 + 



pai^.v) + pn{y,^i) = PD(0,/i) = log' 
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If pn{0, fi) > logM we require po{i^,p) = logM and otherwise we set = 0. In any case 
we always have 

MK^ = ei°g^'^e^°(°''^) <K^ + M. 

It follows that 

(27) [ e^^^^" < e^^' [ e^^'^ < oo, 

so that u satisfies the hypotheses of Theorem 14.11 as soon as we choose M = ^ . We can 
therefore solve the Beltrami equation for u = v{z) to get a principal mapping F of class 
iyi'2(D). Next, set 

2 



f28) 



K[W] 



1- p{z)iy{z) 



w = F{z), z G 



According to Corollary 14.31 n is well defined almost everywhere. We also see that 



1 + 


K, 


1 + 


1 — flU 


1 - 


K 


1 - 





< M 



< oo. 



P 



Thus we may solve the uniformly elliptic Beltrami equation for n to obtain a M-quasiconformal 
principal mapping g. 

We next put f = g o F. The classical Gehring-Lehto result [13J on Sobolev homeomor- 
phisms verifies that F is different iable almost everywhere. Since F has the Lusin property 
jV^^, the same is true for /, and we have 

\Df{z)\^ < MJ{w, g)K{z, F)J{z, F) = MK,{z)J{z, /), w = F{z). 

Arguing as in ([9]) we see that 

\Df\' _ 2 



(29) 



< 



\og{e +\Df\) - p 



is locally integrable. 

In the proof of Theorem 14.11 we constructed the approximants F"^ of the principal 
mapping F. To obtain uniform bounds we may apply the classical Bieberbach area 
theorem [TTj. This gives 



J{z, (j)) < vrr^ 



r > 1, 



B(0,r) 



for any quasiconformal principal solution, conformal outside the unit disc. In particular 
with (|29|1 we see that the compositions goF^ are uniformly bounded in the space Wl^^{<C), 

+2 

Qit) 



log(e + 1) ■ 

Since go F"^{z) f{z) locally uniformly, we infer that / is a Sobolev mapping contained 
in the class Wl^^{c). However, note that although the quasiconformal map g will lie [2] 
in the space W^A'c (C) for all 



2 < s < 



3 — p' 
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yet the composition f = g o F will not lie in Wl^^{c) in general. 

We defined the Beltrami coefficient k, = fig through the formula fl28l) . but as well one 
may identify fig = /i/oF-i via the familiar formula for the dilatation of a composition of 
mappings, 

, , IJ'fiz) — Uf{z) I (z) \2 

1- fif{z)flFiz)^\Fz{z)\/ 

which is a direct consequence of the chain rule. Comparing the expressions shows that 
fif = fi, and hence / in fact solves the required Beltrami equation /j = /i/^. Clearly 
/ is a principal solution. Therefore we have established the existence of the measurable 
Riemann mapping, for any exponentially integrable distortion. 

Uniqueness and Stoilow factorization now follow from Theorem 14.21 Thus the proof of 
Theorem 11.31 is complete. □ 

Once we have the existence of the principal solution, with the equicontinuity properties 
provided by Lemma 12.11 and Corollary 12. 3[ the usual normal family arguments quickly 
give solutions to the global degenerate Beltrami equation, where the coefficient fi{z) is 
not necessarily compactly supported. For an overview, see [3], Section 20.5. 

From the above proof and (1271) we distill a powerful factorization, one of the key facts 
in obtaining the sharp regularity. 

Corollary 4.4. Suppose the distortion function K = K{z) satisfies G L^{IS)) for some 
p > 0. Then for any M > 1 the principal solution to fz{z) = fi{z) fz{z) admits a 
factorization 

where both g and F are principal mappings, g is M -quasiconformal and F satisfies 

Remark. Using the above factorization we obtain the properties M and for all 
maps of exponential distortion. Namely, assume that G W^^^ic) is a mapping of finite 
distortion with e^^^''^'^ G for some p > 0. Then ffist use locally the Stoilow 

factorization (p = ho f where h is holomorphic and / is a principal solution as in Theorem 
11.31 Applying Corollary 14.41 we can make a further factorization, f = g o F where g 
is quasiconformal and F has exponential distortion G with p > 2. According to 
Corollary 14.31 each factor in (p = ho g o F has the properties A/"^^. 

5. Optimal regularity: proof of Theorem 11.11 

We start with general area distortion bounds that have independent interest, since 
they are optimal up to estimates at the borderline case. In fact, the factorization method 
described in Corollary 14.41 enables us to improve the distortion exponent of Corollary! 



Theorem 5.1. Let \fi{z)\ < 1 almost everywhere and fi{z) = outside the unit disk. 
Suppose f is a principal solution to /j = fif^. If 
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then for any < P < p we have 

(30) \f^E)\<C\og-^ie + j^^), EC J}. 

The constant C above depends on (3, p and ||e-'^*^^'-'^^||p only. 
Proof. Choose (3q G and M > 1 so that 

We will then use the factorization f = g o F from Corollary 14.41 Since pM > (3oM > 2, 
Corollary 13.21 applies to cr^ = and F^ = 1 + SF^ = 1 + iScr^, 

\F{E)\= I \F,\'-\F,\'<2\E\ + 2 f |F, - < C log2-*^(e + -^). 

Above is legitimately obtained by integrating the Jacobian since F is a Sobolev 

homeomorphism that satisfies condition A/". On the other hand, since g is an M-quasiconformal 
principal mapping we can use the area distortion estimate f|T6|) . This gives 



\fiE)\ = \goF{E)\ < 7TM\FiE)\'/'' < vrC[log(e + 



(2-/3oA/)/M 



Since P < Pq — the result follows. □ 

Again the family ([3]) shows that the measure distortion bounds are optimal in terms 
of the exponent of the logarithm in (!30|) . up to the possible estimates at the borderline 
P = p. 

We are ready for 

Proof Theorem 11.11 Assume that we have a mapping of finite distortion / G Wl^^{Q) 
with e^^^'f^ G L^^^i^). From it follows that / G VT^^f (f^), and hence the general 
Stoilow factorization Theorem 14.21 applies. With the factorization we may in fact assume 
that Q = C and that / is the principal solution of Theorem II. 3[ conformal outside |d. It 
is then enough to consider the regularity over D. 
We prove first 

(31) J{z, f) log^ (e + J{z, /)) G LL(ro) < /3 < p. 

In order to verify this it is enough to show that J*(x) log^(e + J*{x)) G i^^(0, vr), where 
J*(x), X > 0, is the nonincreasing rearrangement of J(2, /). 
For the rearrangement the estimate fl30l) holds in the form 



(32) 



/■* 1 

/ r{x)dx <C\og-'^{e + 0<p<p, 0<t<n. 
Jo ^ 



Observe that we have applied here the remark after Corollary 14. 4[ Further, since J* is 
nonincreasing, 

(33) rit) < - [ J*{x)dx < -. 

t Jo t 
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Let US now fix < /3 < p and choose a G If we write 

0(t) := f J*{x)dx, < t < vr, 
Jo 

then < (j){t) < Clog""(e + 1/t), 0<t<n. In addition, ([33]) gives 

/ J*{x) log^(e + J*{x)) <C (f)'{x) Iog'^(e + 1/x) dx. 
Jo Jo 

Now an integration by parts shows that the last integral is finite. We have thus proved 

the claim fl3Tl) . 

The second claim 

(34) \Df\'\og^-\e+\Df\)eLlM 



can be deduced from (I3T1) by observing that for every > there are positive constants 
Ci and C2 such that 



xy I 



log^-^(e + y/xy) <CiX log^(e + x) + C2 e^^ for all x,y>0. 



For this consider separately the cases where x < e^^ and where x > e^^. Now, since 
s^s^ log^-i(e + s) is increasing, we conclude 

\Df{z)\Hog^-'{e+\Df{z)\) < K{zJ)J{zJ) log^"^ (e + ^ K{zJ)J{zJ) 



< C,J{z, f) log^ (e + J(^, /)) + C2 eP^(^'^) G LLW- 

Lastly, the family (7^ from ([3]) shows that (|3T1) . (IMl) may fail at the borderline P = p. The 
proof of Theorem 11.11 is complete. □ 

6. Divergent Neumann series 

We next verify that our decay estimate from Theorem 13. H is essentially the best possible. 
One way to do this is simply to observe that if the decay of the L^-norm of the n-th term 
would be of the order 0((n + 1)"^^) with a uniform 6 > 1/2, our proof above would yield 
an area distortion result which would be 'too good', contradicting Example ([3]). We leave 
the details to the reader. Instead, we present here a more concrete approach that is of 
independent interest. 

We start by constructing families of mappings. Suppose first that 7 : (0, 1] ^ C is 
continuous with 

|7(s)|<l for < s < 1 

Then define 

1 + 7(5) ds 



(35) p{t) = p^it) := exp 

or equivalently, 



1 — 7(5) s 
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In any case, we see from ( 135|) that 

t ^ \p{t)\ is strictly increasing on (0,1]. 

If, in addition, 

then p{t) ^ as t ^ 0. In particular, in this case the mapping 

f{z) = ^pi\z\), \z\<l, with fiz) = z, \z\>l, 
\z\ 

defines a homeomorphism of C. 
Next, let a > with 



p(t)=(log^V" that is, 7(t) 



tj ' 'V / a + log(5/t)' 

Consider then the holomorphic family of mappings 



2; 

fxi^) = T-iP\i\z\), Px ■= P(A7), where |A| < 1. 



In other words. 



px{t) = exp 



* 1 + A7(s) ds 



0<t<l, and px{t) =t, t > 1. 



h 1-A7(s) s 

The mappings f\ are all quasiconformal in the entire plane. In fact, they define a holo- 
morphic motion 

$:©xC^C, ^{\z) = fx{z). 

At the boundary, when \ ^ Q E v^e attain well defined mappings fc^ of finite 
distortion. If fact, calculating the complex dilatation shows that \p{z) \ = \'y{t) \ for 1^1 = t. 
Hence the distortion function 

K{zJ^) = \±^, t = |;^| and Id = 1, 
1 - |7(t)| 

which is independent of (, with 

i^T ~ — log — as |z| — * 0. 

a \z\ 

Thus given e > we can choose a < | so that e^^^^ E L^^^(d). On the other hand, one 
computes that 

— 4a 

-^{z) ^ ^ log — 
oz \z\ V \z\ 



-^{z) ^ — log— as \z\ 0. 



Therefore 9^/^ ^ L (c) whenever C belongs to the non-degenerate interval 

II + 
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In particular, A dzf\ does not belong to the Hardy space iJ^(D; L^(c)), and hence the 
sequence of norms || (/x5)"/i||2 is not square summable, i.e. 



and the decay given by Theorem 13.11 cannot be improved at p = 1 . 

7. Applications to Degenerate Elliptic PDE's 

The optimal regularity established in Theorem 11.11 obviously has a number of basic 
consequences, for instance towards removability of singularities. However, our aim here 
is not to consider these consequences systematically. We only indicate one such example, 
namely an application to degenerate elliptic equations. Degenerate elliptic PDE's arise 
naturally in hydrodynamics, nonlinear elasticity, holomorphic dynamics and several other 
related areas. In two dimensions these equations are intimately related to the mappings 
of finite distortion. 

Suppose that u G Wl^liVl) is a (distributional) solution to the equation 

(36) DivA(z)Vu = 0, infiCMl 
Let us assume that u has finite energy, that is 

(37) / (yu{z), A{z)Vu{z)) dm{z) < oo. 

Jn 

In case the domain Q is simply connected, the solution u admits a conjugate function, 
defined by the condition 

(38) Vf = JA{z)Vu. 

Here J{x,y) = {—y,x) is the rotation by a right angle. 
In the complex notation, the function 

f = u + iv 

solves the R-linear equation 

where the coefficients fi{z), v{z) are explicit functions of the elements of A{z), independent 
of the particular solution u. Conversely, for these coefficient /i and v, the real part of any 
solution f = u + iv to (15^ solves the second order equation (15^ . 

Next, let us assume that for almost all z, the matrix A{z) is elliptic and symmetric. 
Then there is a function 1 < K{z) < oo, finite for almost all z, such that 

(40) J^\h\'<{h,Aiz)h)<Kiz)\h\', heu'. 

In terms of the mapping f = u + iv, the conditions (l38l) . ( l40i) take the familiar form 

\Df{z)\'<K{z)J{z,f), 
where the Jacobian can alternatively be written as 

J{z,f) = {Vu{z),A{z)Vu{z)). 



OPTIMAL REGULARITY FOR PLANAR MAPPINGS OF FINITE DISTORTION 



21 



In particular, f = u + iv is a mapping of finite distortion if and only if the solution u 
has finite energy in the sense flTTI) . As an immediate consequence of Corollary 11.21 we find 
optimal conditions guaranteeing the L^-regularity of the derivatives of the solutions u. 

Theorem 7.1. Suppose u G Wl^liVt) is a solution to the equation (!36|) with 

{S/u{z),A{z)Wu{z)) dm{z) < oo. 

// the ellipticity hound K{z) of the coefficient matrix A{z) satisfies 

(41) e^(^) G LL(fi) 
for some p > 1, then the solution u has the regularity 

(42) \Vu\'' + \A{z)Vu)\^eLl,{n). 

Similarly, we have 

I Vm|2 \ogf'-\e + \Vu\) G Li^iSl), 0<(3<p, 

and 

\A{z)Vu)\' \og^-\e + \A{z)Vu)\) G LUQ), < f3 < p 

whenever the ellipticity K(z) satisfies e^^^-* G Lf^^(^l), p > 0. Furthermore, via the Stoilow 
factorization many properties of the harmonic functions generalize to the solutions to these 
degenerate elliptic equations. There are even non-linear counterparts, see [1]. 

One of the fundamental features of the mappings of bounded distortion are their self- 
improving regularity properties. Similar phenomena occur also in the theory of mappings 
of finite distortion. A typical example is obtained by combining the Stoilow factorization 
Theorem 14.21 with our Theorem II. 1[ There are factorization theorems slightly beyond 
the class W^^^ic), Q{t) = t'^\og~^{e + t), see for instance [H], p. 280. Both factoriza- 
tion methods lead to removability results for mappings of finite distortion. We refer the 
reader to [13], [12], [3] and their references for these and other related methods, such as 
Caccioppoli-type estimates. In view of the discussion in this section, one obtains corre- 
sponding removability results for the solutions of the elliptic PDE's. 

Acknowledgements. The authors want to thank Albert Clop for pointing out a mistake 
in an earlier version of Example 
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